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Abstract 

A modified method for calculating the non-perturbative quark vacuum condensates from the 
global color symmetry model is derived. Within this approach it is shown that the vacuum con- 
densates are free of ultraviolet divergence which is different from the previous studies. As a special 
case the two quark condensate (qq) and the mixed quark gluon condensate giqG^^o'^'^q) are cal- 
culated. A comparison with the results of other nonperturbative QCD approaches is given. 
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I. INTRODUCTION 



The non-perturbative structure of the QCD vacuum is characterized by various con- 
densates: such as the two quark condensate (gg), the mixed quark-gluon condensate 
g{qGi^,^a'^'^q) , the four quark condensate {qTqqFq), etc. These condensates are essential 
for describing the physics of the strong inter action [1,2]. There has been growing interest in 
describing the properties of hadrons in nuclear matter in terms of the in-medium quark and 
gluon condensates, which are shifted from their vacuum values [3-8]. Since the Global Color 
Symmetry Model(GCM)[9] provides a nonperturbative framework which admits a simulta- 
neous study of the spontaneous chiral symmetry breaking and confinement, it is useful in 
exploring the phase transition at nonzero temperature and density. Before we generalize 
GCM from zero temperature and density to nonzero temperature and density, it is neces- 
sary to check whether GCM can provide a good description of vacuum properties at zero 
temperature and density. 

It is the aim of this paper to consider quark vacuum condensates in general and in particu- 
lar the two quark condensate and the mixed quark-gluon condensate in the framework of the 
global color symmetry model. Although GCM violates local color SU{3)c gauge invariance 
and renormalizability, it provides a very successful description of various nonperturbative as- 
pects of strong interaction physics and hadronic phenomena at low energies. These include, 
for instance, quark confinement [10], ?7a(1) breaking, and the rj — rj' mass splitting[ll], low 
energy chiral dynamics of Goldstone bosons(7r, iC, 77) [12-18], meson form factors[19], heavy- 
light mesons[20], systems at finite temperate [21], or soliton[22] and Fadeev[23] descriptions of 
the nucleon. The first estimations of the two quark and the mixed quark-gluon condensates 
in GCM was made in Ref.[24], where the two quark condensate (qq) itself is systematically 
smaller than the "standard" value of -(250Mey)^[l] 

Here, we present a modified method for calculating the quark vacuum condensate in 
the framework of GCM. Since the u, d current quark mass is small, we consider the GCM 
generating functional for massless quarks, i.e., the chiral limit, in Euclidean space 




(1) 



where 



Sgcm — 



d'^x q{x) 7 • 5 — ig-^'y ■ A"'{x) q{x) + 



d'xd'y ^Alix) [d^^{x - y)Y' Al 
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with D'^'l{x — y) denoting the gluon two-point green function. Because the form of — y) 

in the infrared region is unknown, one often treats the — v) ^ the GCM input 

parameter, which, as we will discuss later, is chosen to reproduce certain aspects of low 
energy hadronic properties. For convenience we will use the Feynman like gauge D'^Hx—y) = 



5fj_y5 D[x — y) from now on. 

Performing the functional integral over DA in Eq.(l)(see Eq.(27) below), we obtain the 
GCM generating functional 



(2) 



A 

here denotes the color octet vector current = <i{x)'-)^^q{x). 

Introducing an auxiliary bilocal field B^{x,y) as in Ref.[9], the generating functional of 
GCM can be written as 



where 



Z[ri,'n] = j '^Q^Q'^B%x,y)exp^-S[q,q,B%x,y)]+ j d'^x{r]q + qr})^ , 
S%q,B\x,y)]= j jd'xdS \q{x)g-'[x,y- [B^My) + ^'/"".'^jf l^^^"^^ 



(3) 



with 



1 



g-\x, y- [B%] = 7 • dj('^ {x-y) + -A'B\x, y) , 



(4) 



where the matrices = D"'F^C'^ is determined by Fierz transformation in Dirac, flavor and 
color space, and are given by 



1 1 



A ± 
^3 ^'^3 



Performing the functional integral over "Dg and Vq in Eq.(3), we obtain the GCM gener- 
ating functional 

Z% ri]^l VB%x, y) exp(-5[77, 77, B^x, y)]) , (5) 

where 



S\r]^ri^ B^ {x^y)] = — Trln 

+ J J d^xd^y 



p5{x-y) + -k'B\x,y) 
B\x,y)B\y,x) 



2gW{x-y) 



+ fj{x)g{x,y;[B'My) 



(6) 
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The saddle-point of the action is defined as 5S[r], rj, B^{x, y)]/5B^{x, y) =0 and is 

77=77=0 

given by 

Bl{x -y) = g^D{x - y)tr[A'go{x - y)], (7) 

where Go stands for ^[-Bq] and the trace in Eq.(7) is to be taken in Dirac and color space, 
whereas the fiavor trace has been separated out. 

We will calculate the vacuum condensates from the saddle-point expansion, that is , we 
will work at the mean field level. This is consistent with the large Nc limit in the quark 

fields for a given model gluon two-point function. In the mean field approximation, the field 
B^{x — y) is substituted by their vacuum B^^x — y). Under this approximation, the dressed 
quark propagator G{x — y) = Qq{x — y) in GCM has the decomposition 

G-\p) = i^-p + T.{p) = i7 • pAip") + Bip") (8) 

with the self-energy dressing of the quarks E(p) is defined as: 

E(p) ^ \h'Bl{p) = J d'xe^^-^ [l^'B'.ix)] = n ■ p[A(p') - 1] + B(p% (9) 

where the self energy functions A{p'^) and B{p^) are determined by the rainbow Dyson- 
Schwinger equation [9] 



f) 

In order to get the numerical solution of AijP') and B^p"^), one often uses model forms for 
the gluon two-point function as input in Eqs. (10-11). Here we investigate two different two 
parameter models for gluon propagator 

g^D(^\q^) = g-D^l{<f) + g^DMl') = 37r^|^e-^ + ^ ^, (12) 

\'^QCD J 

and 

g-D^\e) = g-Dfli^) + g^Duv{q') = ^^'d-^ + /^'^ . . (13) 

The term DiR{q'^), which dominates for small q^, simulates the infrared enhancement 
and confinement. The other term Duvicf) dominating for large (f is an asymptotic ul- 
traviolet (UV) tail which matches the known one- loop renormalization group result with 



d = [12/(33 - 2A^/)] = 12/27, Aqcd = 200 MeV. The model parameters x and A are 
adjusted to reproduce the weak decay constant in the chiral hmit /tt = 87 MeV. The forms 
of g'^D{q^) have been used in Ref.[15] and it has been shown that with these values a satis- 
factory description of all low energy chiral observables can be achieved(more detail can be 
seen in Refs.[15] and [24]). 

Here we want to stress that the B{p^) in Eqs.(10,ll) has two qualitatively distinct solu- 
tions. The "Nambu-Goldstone" solution, for which 

^ 0, (14) 

describes a phase in which: 1) chiral symmetry is dynamically broken. Because one has a 
nonzero quark mass function; and 2) the dressed quarks are confined, because the propaga- 
tor described by these functions docs not have a Lehmann representation [10]. In "Nambu- 
Goldstone" phase, the vacuum configuration Bq[x — y) in GCM(at the mean field approx- 
imation) can be regarded as a good approximation to the "exact" vacuum in QCD. The 
alternative "Wigner" solution, for which 

B{p^) ^ 0, (15) 

describes a phase in which chiral symmetry is not broken and the dressed-quarks are not 
confined. In "Wigner" phase, the vacuum configuration Bq{x — y) in GCM(at the mean 
field approximation) corresponds to the "pcrturbative" vacuum in QCD. 

With these two "phase" characterized by qualitatively different quark propagators, the 
GCM can be used to calculate the vacuum condensates as we will show later. 

Our paper is organized as follows: In Sec. II we briefly review the definition of vacuum 
condensate in QCD sum rule, and exhibit a modified method, which is consistent with this 
definition, for calculating the non-perturbative quark vacuum condensates from the GCM. 
Based on the modified approach, the two quark condensate and the mixed quark gluon 
condensate are calculated at the mean field level. A brief discussion and conclusion is given 
in Sec.III. 

II. EVALUATION OF VACUUM PROPERTIES 

In order to ensure that the treatment about the vacuum condensates in this paper is 
consistent with that in QCD sum rule, a brief introduction of vacuum condensates in QCD 
sum rule is described below: 
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In QCD sum rule, one often postulates that quark propagators are modified by the long- 
range confinement part of QCD; but the modification is soft in the sense that at short 
distance the difference between exact and perturbative propagators vanishes. 

To formalize this statement, one can write the "exact" propagator G{x) as a vacuum 
expectation of a T-product of fields in the "exact" vacuum |0): 



G,,{x,y) = {0\T[q,{x)q,{ym. 
According to the Wick theorem, one can write the T-product as the sum 

T[qi{x)qj{y)] = qi{x)qj{y)+ : qi{x)qj{y) : 



(16) 



(17) 



of the "pairing" and the "normal" product. The "pairing" is just the expectation value of 
the T-product over the perturbative vacuum |0) 



q,{x)qj{y) = (0|Tfe(x)^,(y)]|0) ^ G',f\x,y), 



(18) 



i.e., the perturbative propagator. By this definition, we have the following two quark vacuum 
condensate (0| : q{x)q{y) : |0): 



(0| : qi{x)qj{y) : |0) = {0\T[q,{x)qj{ym - {0\T[q,{x)q,{y)]\0) 

= (-) \Gj,(y,x) - G',r(y,x)] ^ H^Ay,^) = (-) 



(19) 



G{q^) - GP'^-'iq^) 



where E(a;, y) — G{x, y) — G^^''*(a;, y). Thus, our assumption of (0| : gg : |0) 7^ is equivalent 
to the statement G{x) ^ G^^''*(x). Because at large momentum region the difference between 
exact and perturbative quark propagator vanishes, the quark condensate (0| : gg : |0) is free 
of UV divergence. This conclusion is reasonable because the vacuum condensates reflects 
the IR behavior of QCD. 

Similarly, we have the four quark vacuum condensate: 

(0| : g(a;)A(^)g(a;)g(t/)A(2)g(t/) : |0) 

= (0|r [q{x)K^^^q{x)q{y)K^''\{y)] |0) - (0|T [q{x)K^^^ q{x)q{y)K^^\{y)] |0) 
-(0| : g(x)AWg(x)g(y)A(^)g(y) : |0) - (0| : qix^h'''^ q{x)q{y)A'^ : |0) 



(20) 



-(0| : g(x)AWg(x)g(y)A(^)g(y) : |0) - (0| : q{x)h.^^\{x) q{y)A.^''\{y) : |0), 



here the A^'^ stands for an operator in Dirac and color space. It should be noted that 
the treatment of vacuum condensates here is different from that in Refs. [24,25], where the 



contribution of the first term of the right hand of Eq.(20) is considered to be the four quark 
vacuum condensate. By this definition, the vacuum condensates in Refs. [24,25] are just the 
corresponding quark green function taken at one point. This is not consistent with the 
definition of quark vacuum condensate in QCD sum rule. 
In addition, we have the six quark vacuum condensate: 

(0| : q{x)A^^\{x)q{y)A^^\{y)q{z)A^^\{z) : |0) 

= tr[E{y, y)A^^'^]tr[E{z, x)A^^'>E{x, z)A^^^] - tr[E{z, x)A^^^E{x, y)A^^^E{y, z)A^^^] 
-tr[^{y, x)A«E(x, z)A^^^^{x, y)A^^^] + tr^x, x)A«]tr[E(x, y)A^^^^y, z)A^^^] (21) 

+tr[E{z, z)A^^'>]tr[J:{y, x)A^^^J:{x, y)A^^^] - tr[S(a;, x)A^^^]tr[J:{y, y)A'^^^]tr[J:{z, z)A^% 

Based on the above statement, in order to calculate quark vacuum condensates one must 
know not only the "exact" but also the "perturbativc" quark propagator in advance. The 
calculation of the "exact" quark propagator in GCM(at the mean field approximation) has 
been given by Eqs. (10,11) and Eq.(14), the left question now is how to treat consistently 
the "perturbative" quark propagator in GCM. 

In "Wigner" phase(the perturbative quark scalar self energy function B'{p^) = 0), the 
Dyson- Schwinger equation (10, 11) reduces to: 

(22) 

where A^p^) denotes the perturbative quark vector self energy function. Therefore, the 
perturbative quark propagator in GCM can be written as GP^^{q'^) — j^^q2)q2 — —il-(lC{q^). 
Numerical studies show that for q^ » -^qcd) has 

oi^) - cr-Xe) - -'2^^l^ + n ■ ,c(,^) - 0, (23) 

which is just what one would expect in advance. The above results can be seen from the 
following figures (fig.l-fig.4). 
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Fig.l. Comparison between CJ^(s) and C(s) 



Fig.2. Comparison between CJ^(s) and C(s) 



- (5^(s) set 1 (A=0.2GeV^, 3C=l-65GeV) 
C(s)setl 

<5^(s) set 2 (A=0.02GeV^ x=l-55GeV) 
-- C(s)set2 

- - <5^(s)set3(A=0.002GeV^,x=1.45GeV) 

- - C(s)set3 




C(s) set 1 

<5^(s) set 2 (A=10"'GeV'', x=l-02GeV) 
CCs) set 2 



C(s)set3 

c.(s)=A(s)/(a'(s)s+b'(s)) 




1 1.5 

s (GeV^) 



1 1.5 

s (GeV^) 



2„(1) 



Fig.3. 0„(s) for gluon propagator g D 



2„(2) 



Fig.4. 0„(s) for gluon propagator g D 



0.5 - . 



set 1 (A=0.2GeV , x=1.65GeV) 
set 2 (A=0.02GeV^, x=1.55GeV) 
set 3 (A=0.02GeV^, x=1.45GeV) 

ct„(s)=B(s)/(a'(s)s+B%)) 




~l ' \ ' 1 ' — 

set 1 (A=10"'GeV^ x=l-83GeV) 

- set2(A=10"'GeV'', x=1.02GeV) 
set 3 (A=10"'GeV'', x=0.83GeV) 

CTg(s)=B(s)/(A'(s)s+B%)) 




Once the "exact" and "perturbative" quark propagator in our model are determined, one 
can calculate the (0| : qq : |0), (0| : qK^^^qh.^^^ q : |0) vacuum condensate and the mixed 
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quark gluon condensate g{qGfj,i,cr^^q) at the mean field level. In particular we obtain the 
two quark condensate (0| : qq : |0) in the chiral limit; 



(0| : : |0) = (-) lim ir[E(x,0)] 

(-) lim tr[G(x, 0)] = (-) lim / (2^),^ ^.^.^^.^ ^ ^^^.^ 
(-) / JPl^i^(^ = A r rf, ^^(^) (24) 



where we have used the fact that tr[GP^^*{x, 0)] = in the chiral limit. Due to the contribu- 
tion of the subtraction term vanishes, this is the quark Green function taken at one point, 
which in momentum space is a closed quark loop. Therefore, in principle, the upper limit 
of integration in Eq.(24) should be taken to be infinity. 

It should be noted that the upper limit of integration in Eq.(24) is different from that in 
Refs. [24,25]. In Refs. [24,25], the two quark condensate is calculated as 

where fi is interpreted as the renormalization scale in Ref.[24]. This is incorrect (/x is only a 
hard cutoff) . Due to the misuse of the cut off in the above integration (/x is chosen to be 
IGeV'^ in Ref.[24]), the value of two quark vacuum condensate is strongly underestimated in 
Ref . [24] (see Table III and IV below) . Here we stress the fact that one found from perturbative 
theory[26] 

+ 



B{j?) = mo ( 1 — ^^In 



p2 
ml 



47r 

i.e., the perturbative correction is proportional to the current quark mass m^. Hence in a 
perturbative analysis B{p'^) oc mo, which vanishes in the chiral limit mo 0. This means 
that the nonzero B{p'^) even in the region of p^ > IGeV^ is an intrinsically nonperturbative 
effect in the chiral limit case. Therefore, the following contribution 

3 , sB(s) 

as- 



=iGey2 sA'^{s) + B{s) 

should be taken into account in a correct treatment of the two quark vacuum condensate, 
that is, the upper limit of integration /j, should be taken to be infinity. 

Another important vacuum condensate is the four quark condensate which are needed 
within the framework of QCD sum rule to describe properties of both mesons and baryons. 
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Prom Eq.(20), it is easy to obtain 

(0| : q{x)M'^q{x)q{y)M^^q{y) : |0) 

= - {ir[E(y,x)A«E(x,y)A(^)] - ir[E(x, x)A«]ir[E(y, y)A(2)]} . 



(25) 



By means of Eq.(25), one can calculate all kinds of four quark condensates at the mean field 
level in GCM. For instance, in the case of A*^^) = A^^^ = 75 and x — y — 0, one finds 



(0| : q{Oh,q{0) g(0)75g(0) : |0) 

d^p d^q \\ B{p^] 



(-12)/ 



+ 



(27r4) (27r4) 
A{p') 



A^{p^)p^ + B^{p^) 
A{q' 



B{q' 



A^{p^)p'^ + B^{p^ 



A^{q^)q^ + B\q^) 



A\q^)q^ + B\q' 



p-q. 



(26) 



Notice, the integral in Eqs.(26) is free of UV divergence. However, if we follow the method 
proposed in Ref.[24], we have 



: ^-(0)759(0) g(0)75g(0) : \0)^~ 
d^p A [[ B{p^) 



(-12)/ 



+ 



(27r4) (27r4) 
A{p') 



A-^{p^)p^ + B^{p^] 



A{q') 



B{q' 



A\q^)q^ + B\q^) 



_A-'{q^)q^ + B\q^) 

which is UV divergent. Here it should be noted that the calculation of vacuum condensates 
in the framework of GCM is free of UV divergence if one adopt the adequate subtraction 
mechanism. This result is also quite different from that of effective quark-quark interaction 
model, such as Nambu, Jona-Lasinio[5] and chiral constituent quark models[27]. 

The mixed quark gluon condensate g{qGfj,u(j'^^q) is also an important condensate within 
the framework of QCD sum rule. Now let us focus on the gluonic obscrvables. Because the 
functional integration over the gluon field A in Eq.(l) is quadratic for a given quark-quark 
interaction D we can perform the integration over gluon fields analytically. Using the short 
hand notation for the integrations as in Ref.[24], we have 



/ 



VA e-5^^-^^+M = e^^'^^ I VA ^(x)e-5^^-'^+^-^ 



(jD)(x)e^^-^^ 



j VA A{x)A{y)e~^ 



AD-^A+jA 



D{x,y) + {jDf{x,y) 



62 



(27) 



This means that the gluon vacuum average renders effectively a quark color current 
q{x)^fj,-§-q{x) together with the gluon two-point function D. In order to calculate the mixed 
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quark gluon condensate g{0\ : qG f^ucr'^^ q |0), we have to do integration over powers of gluon 
field A and A^. By means of Eq.(27) we find 



^(0| : q{x)G,,{x)a'^''q{x) : |0) 

2g{0\ : qix)d;Atix)^a^''q{x) : |0) +/r^^(0| : qix)A'^ix)Atixy-^a^''q{x) : |0) 



= 2ig^ J d^z [d^^^D{x,z) 



■.q{x)q{z)^^^q{z)^a>^'^q{x):\Q) 



(28) 



A" 



- r^^y d'z,d'z2g'Dix,z^)g'Dix,Z2){0\ : g(x)g(^i)7.yg(^i)g(^2)7.y9(^2)ya^^g(a;) : |0), 

with /"^"^ the structure constants of the SU(3) color group. As indicated by Eq.(28), in order 
to calculate the mixed quark gluon condensate, one must calculate the four quark and six 
quark condensates in advance. Applying Eqs.(20) and (21), we have 



g{0\ : q{x)G^,(x)a^''q{x) : |0) = -8ig' J d^z 

+12i j d'^zid'^Z2g'^D{x,zi)g^D{x,Z2)trD [S(z2, a;)a^j.a(a;, 2;2)7^S(zi, ^2)7^ 

d'^p B{p^)[2 - A{p^)]p^ 81 f d^p B{p^)[B\p^) + 2p^A{p^)-l)A{p^)] 

4 



= -36 



(27r)4 A2(p2)p2 + ^2(p2 



(27 

9 r d^p 
^2 J (2^ 



(2vr; 



16 - 9A{p 



A^{p^)p^ + B^{p^) 



[9A(p2) 



A^{p^)p^ + B'^ip^] 
B(p')C(pV, 



+ 9C(/) 



Bip')p . qCiq') 



(29) 



where tr^ denotes that the trace in Eq.(29) is to be taken in Dirac space only, whereas the 
color trace has been separated out. It is easy to find that our result for the mixed quark 
gluon condensate is substantially different from that in Ref.[24](see Eq.(25) in Ref.[24]). The 
result for the mixed quark gluon condensate in Ref.[24] keeps only the contribution of the 
first and the second term of the right hand in Eq.(29). 
By means of Eq.(22), Eq.(29) can be further reduced to: 



^(0| : q{x)G,^{x)a^^q{x) : |0) = | / ^,B{p') 



^/C(/)+^(P>W)-^ 



9 

+- 



d'^p 



B{p^ 



18A(p2 



32 



4 J (27r)4/12(p2)p2 + 52(p2 

81 r°° 

/ sdsB(s) 
Jo 



327r2 
9 



C(p2 



-sC{s) + A{s)sC{s) - - 



A(p')p' 9A(p') - 16 



647r2 



B(s) ( 18A(s) -32 , r X 



(30) 
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TABLE I: The two quark and mixed quark gluon condensate for g'^ Dj^{q'^)=3TT^ ^ 

AjGeV^] x[GeV] (-)(0| : qq : |0)^ (MeV) (-) ff (gG^.cr^'-g) 5 (MeV) 

0.200 1.65 252 296 

0.020 1.55 212 353 

0.002 1.45 190 349 



TABLE II: The two quark and mixed quark gluon condensate for g'^DjJi^{q'^)=ATr'^d-^ 
A[GeV*] x[GeV] (-)(0| : qq : |6)i (MeV) (-) g{qG^,af'''q)-^ (MeV) 



10-1 


1.83 


344 


301 


10-^ 


1.02 


264 


336 


10-7 


0.83 


264 


393 



Prom Eqs.(24) and (30), we can calculate the two quark condensate and the mixed quark 
gluon condensate. In Table 1 and 11, we display the numerical results of the two quark and 
the mixed quark gluon condensate for three different sets of parameters of two different 
model gluon propagators (Eq.( 12) and Eq.(13)). 

Table I and II show that the result of the two quark condensate and the mixed quark 
gluon condensate is compatible with the corresponding "standard" value determined from 
QCD sum rule [1,2]. Here we want to stress that the form of gluon propagator and the 
corresponding parameters are somewhat flexible, as indicated by the authors in Ref.[28]. 
Thus our result about the two quark and the mixed quark gluon condensate may provide a 
further constraint on the the form and parameters of the gluon propagator. 

It should be noted that the cut-off fj, plays an essential role in calculating the value of the 
condensates in Refs. [24,25]. However, as is shown above, in a correct treatment of vacuum 
condensate, the cut off should be taken to be infinity. This situation is very similar to the 
determination of vacuum condensates in the instanton liquid model where the upper limit of 
integration /j, is taken to be inflnity[29,30]. Because the calculation is numerical, one has to 
use a very large but finite value of upper limit of integration. In our numerical calculation, 
the upper limit is chosen to be 500 GeV'^, which is large enough to ensure that the calculated 
value is independent of the choice of the upper limit of integration. This result can be seen 
from Table HI and IV. 
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TABLE III: The sensitivity of the vacuum condensates to the choice of the upper hmit of integration 
fi in Eqs.(24) and (30) for g'^Dfl{q^)=3Tr^^ei {A = 0.020GeF^, x = 1.55GeF) 



/i {GeV^) (-)(0| : qq : |0)s {MeV) {-)g{qG^,,a>'-q)\ {MeV) 
1 170 352 
100 205 353 
300 210 353 
450 212 353 
500 212 353 



TABLE IV: The sensitivity of the vacuum condensates to the choice of the upper hmit of integration 
II in Eqs.(24) and (30) for Dfl{q^)=ATr'^ d^^{/\ = m'^GeV^, x = 1.02Gel^) 

/X {GeV) (-)(0| : qq : |6)5 {MeV) (-) g{qG^,a^'''q)l {MeV) 



1 177 333 

100 254 336 

300 262 336 

450 264 336 

500 264 336 



III. DISCUSSION 

To summarize, in the present paper, we provided a general recipe to calculate the vacuum 
condensates at the mean field level in the framework of GCM. This approach is quite different 
from that in previous studies[24,25] in the following aspects: 

1. The first difference between the approach proposed in Refs. [24,25] and our's comes 
from the definition of vacuum condensate. The definition of the vacuum condensate in 
Rcfs. [24,25] is only the corresponding Green function taken at one point in coordinate 
space(see Eq.(19) of Ref.[24])), which is not the expectation value of the "normal" 
product over the vacuum |0). Therefore, it is not consistent with that in the QCD 
sum rule approach and the spurious contribution of the perturbative terms have not 
been subtracted. In the calculation of the two quark condensate, the contribution 
of the subtraction term vanishes which explains why our formula for the two quark 
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condensate exactly coincides with the corresponding ones in Ref.[24](in this case, the 
main numerical difference between the two parties comes from the integration region 
between fi—l GeV^ to infinity). In other cases, the contribution of the corresponding 
subtraction terms do not vanish, and this will cause the substantial difference between 
the expression of the vacuum condensates in Ref. [24] and in ours(for example, see the 
expression of the mixed quark gluon condensate in Eq.(29)). In the present paper, 
an adequate subtraction mechanism which is consistent with the definition of vac- 
uum condensate(the vacuum expectation value of the "normal" product of fields) in 
QCD sum rule is adopted to calculate vacuum condensates (employing Wick theorem 
and considering in addition the usual Nambu-Goldstone phase a perturbative Wigner 
phase of GCM). Within this approach all kinds of vacuum condensates arc free of UV 
divergence. This is reasonable and what is to be expected. In addition, due to the fact 
that the definition of the vacuum condensate in Ref. [24] is not the vacuum expection 
value of the normal product of fields, it will be plagued with UV divergence in some 
cases(see Eq.(26) and related items). 

2. The second disagreement between the two parties is the value of the upper limit 
of integration In the calculation of vacuum condensates of the present paper, all 
upper limit of integration is taken to be oo. This is very similar to the determination 
of vacuum condensates in the instanton liquid model [29,30], which is different from 
that in Refs.[24] and [25]. Because of the misuse of cut oS /i — IGeV'^ in Ref. [24], the 
value of the two quark condensate is strongly underestimated there. In addition, if we 
use Eq.(25) in Ref. [24] and let the cutoff // to be oo to calculate the mixed quark gluon 
condensate, we get a value which is several times larger than the "standard" value of 
the mixed quark gluon condensate obtained by QCD sum rule. 

The numerical results of the two quark condensate (qq) and the mixed quark gluon 
condensate g {qG ^vCJ^^'^ q) obtained by the modified approach show that they arc compatible 
with the values obtained from QCD sum rule. These results can be used as a consistency 
check of the GCM and provide a reference for further study. Finally, we want to stress that 
the GCM is not renormahzable. Therefore, the scale at which a condensate is defined in 
our approach is a typical hadronic scale, which is implicitly determined by the model-gluon 
propagator g'^D{q^) and the solution of the rainbow DS equation (10-11). This situation 
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is very similar to the determination of vacuum condensates in the instanton hquid model 
where the scale is set by the inverse instanton size[29,30]. 
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